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Abstract. By fitting the available photometric data on V838 Mon with standard supergiant spectra we have 
derived principal stellar parameters, i.e. effective temperature, radius and luminosity, and followed the evolution 
of the object since its discovery in early January 2002. Our analysis shows that the 2002 outburst of V838 Mon 
consisted of two major phases: pre-eruption, which was observed in January 2002 and a major outburst, the 
eruption, which started in the beginning of February 2002. During pre-eruption the object seemed to be relaxing 
after an initial event which had presumably taken place in last days of December 2001. The eruption phase, which 
lasted untill mid- April 2002, resulted from a very strong energy burst, which presumably took place in last days 
of January at the base of the stellar envelope inflated in pre-eruption. The burst produced an energy wave, which 
was observed as a strong luminosity flash in the beginning of February, followed by a strong mass outflow in 
the form of two shells, which was observed as an expanding photosphere in later epochs. In mid-April, when the 
outflow became optically transparent and most of its energy radiated away, the object entered the decline phase 
during which V838 Mon was evolving along the Hayashi track. This we interpret as evidence that the main energy 
source during decline was gravitational contraction of the object envelope inflated in eruption. Late in 2002, dust 
formation started in the expanding shells which gave rise to a strong infrared excess observed in 2003. 
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1. Introduction 

The eruption of V838 Mon was discovered in the beginning 
of January 2002 (Brown 2002). It appeared to be rather 
a pre-eruption brightening as the main eruption started 
at the beginning of February 2002 and lasted for about 
two months (e.g. Munari et al. 2002c). The rise to the 
main maximum in February 2002 caused to a light echo 
(Henden et al. I2002fl whose spectacular images were ob- 
tained with HST (Bond et al. l2TjO"3|l . As shown in Tylenda 
(2004) an analysis of the echo images allows us to study 
the details of the dust distribution near V838 Mon, as well 
as to constrain the distance to the object. 

Early classifications of the eruption of V838 Mon 
as nova-like quickly appeared to be inappriopriate (e.g. 
Munari et al. I2002cl Soker & Tylenda EES)) . The light 
curve of V838 Mon was atypical for novae. Also, the ob- 
served mass loss seemed to be not enough violent for a 
nova. But the strongest argument against nova-like char- 
acter came from the spectral evolution of V838 Mon. 
Classical novae usually reach a minimum effective tem- 
perature near their optical maximum brightness and then 
tend to evolve toward higher effective temperature, reach- 
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ing values well above 10 5 K. V838 Mon, on the contrary, 
reached the maximum effective temperature, characteris- 
tic of an A-F spectral type, at the optical maximum in 
the beginning of February. The object then showed a gen- 
eral tendency to evolve to lower effective temperatures. In 
mid- April 2002 it almost disappeared from optical images 
but remained very bright in infrared, becoming one of the 
coolest M-type supergiants yet observed in astrophysics. 
Detailed descriptions of the spectral and photometric evo- 
lution of V838 Mon can be found in a number of papers in- 
cluding Munari et al. (2002c), Kimeswenger et al. lP002), 
Kolev et al. (|2002(1 . Osiwala et al. H2003|l . Wisniewski et al. 
(HUH), Crause et al. HUU^i . Kipper et al. l|2Tjtj4"|> . 

The global fading of V838 Mon in optical after eruption 
enabled us to discover a faint hot continuum in short wave- 
lengths (Desidera & Munari 2002 Wagner & Starrfield 
2002 ) later classified as a normal B3 V star (Munari et al. 
2002a} • A detailed analysis of observational data available 
for the progenitor and the enviroment of V838 Mon has 
been carried out in Tylenda, Soker & Szczerba (2005b). 
The main conclusion of this study is that V838 Mon is 
most likely a young binary system consisting of two inter- 
mediate mass stars. One component is the B-type compan- 
ion discovered by Munari et al. (|2002a|l . The other com- 
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ponent that erupted in the 2002 event has a mass between 
5-10 Mq and before eruption was either a main sequence 
star of similar mass and spectral type as the B-type com- 
panion or a slightly less massive pre-main-sequence star. 
(However there is no direct evidence that both compo- 
nents form a real binary.) The system is embedded in the 
interstellar medium, which is most likely to be related to 
molecular clouds seen around the position of V838 Mon, 
rather than to be expelled by the object itself as argued 
by van Loon et al. i|2004|l . 

The aim of this paper is to study the evolution of the 
effective temperature and luminosity of V838 Mon dur- 
ing and after its eruption. This sort of study is crucial 
for testing possible eruption mechanisms and can be an 
important step towards a final understanding of the na- 
ture of V838 Mon. Our analysis is based on photometric 
results published in the literature. This is the only kind 
of observational data for V838 Mon where the available 
results are sufficiently rich to undertake a uniform and 
meaningful analysis over a time span since the discovery 
of the eruption untill the present (December 2004) state 
of the object. 

2. Observational data and analysis 

The analysis done in the present study is based on a com- 
parison of the observed photometric magnitudes to a set of 
reference spectra. The least squares method has been used 
to get the best fit between the model spectrum and the 
observations. Since V838 Mon is considerably reddened 
(see Sect. El) the model spectra have been reddened before 
fitting them to the observations. The standard extinction 
curve, R = Av/Eb-v = 3.1, has been adopted. From the 
best fit one obtains an effective temperature, T e ff, and an 
angular stellar radius, 9. If the distance to the object is 
known, a linear radius, R, and a stellar luminosity, L, can 
then be calculated. 

It is clear that the larger the wavelength range that is 
covered by the photometric measurements, the more reli- 
able are the results obtained. Therefore we have primarily 
based our analysis on UBVR C I C JHKL magnitudes avail- 
able from Munari et al. I|2002c|) and Crause et al. I|2()()3I 
2005}. For certain time periods no JHKL measurements 
are available and we had to rely on UBVR C I C data from 
the above references, supplemented by the BVRI C mea- 
surements of Kimeswenger et al. |2nQ2J). In later epochs 
when the object became very cool, the infrared data were 
particularly important. Therefore we have also used pho- 
tometric data from Henden & Munari (2002), Watson & 
Costero $M)ty . Lynch et al. lj2lK73l l2l?0ljl and Tapia & 
Persi (EE!!- 

Given the luminosity range obtained in the present 
study (10 4 — 10 6 Lq) as well as the spectral appearance 
of V838 Mon during eruption (e.g. Kipper et al. [2004) 
we have used standard intrinsic colours for supergiants 
(luminosty class I) when constructing the set of refer- 
ence spectra. The intrinsic photometric colours have been 
taken from Schmidt-Kaler lfT38"2l - UBV), Johnson 0966 



- VRI for FG K type s), Lee (P7UI - VRI for M types) 
and Koornneef (|1983l — V JHKL). A calibration of the ef- 
fective temperature and bolometric correction against the 
spectral type has been taken from Schmidt-Kaler (1982). 

In this way we obtained a grid of 21 reference spectra 
ranging from F0 1 to M5 1 (T c ff ranging between 7700 and 
2800 K). Between the grid points a linear interpolation 
was done either in the colour - log T e g plane (e.g. B — V 
or V — K versus log T ff) to get a reference spectrum for 
any given value of T e g or in the colour - log A plane if the 
effective wavelength of a photometric band (e.g. R c , I c ) 
differed from the Johonson standards. 

After mid- April 2002 V838 Mon became cooler than 
the M5I standard. In order to have estimates of the ob- 
ject parameters in later epochs we proceeded in two ways. 
First, we extrapolated our grid of standard spectra down 
to M10. The extrapolation is in the log colour - log T e s 
plane (including log BC - log T c g). In the range of M 
types we explore the Wien part of the spectrum and conse- 
quently the relations are roughly linear in the above plane. 
Second, we fitted a blackbody spectrum to the observa- 
tions. For the effective wavelengths of particular photom- 
teric bands, the fluxes obtained from the blackbody law 
have been converted to magnitudes using the standard 
calibration of Vega rSchmidt-Kaler lT5%2l Tokunaga|m3|) . 
After having reddened them according to the extinction 
curve, the magnitudes were compared to the observations. 

3. Distance and interstellar extinction 

Initial estimates of the distance to V838 Mon done in 
Munari et al. ( 2002c) and Kimeswenger et al. 2002) gave 
values of 0.6-0.8 kpc. They were however based on a 
wrong interpretation of the observed expansion of the light 
echo. From a more realistic study of the echo structure 
Bond et al. ( 2003) concluded that the distance is at least 
6 kpc. A detailed study of the light echo evolution done in 
Tylenda PMjl and Tylenda et al. ( 2005b ) shows that the 
distance cannot be unambiguously determined from the 
observed expansion of the outer echo edge. Only a lower 
limit of ~ 5 kpc can be established. However, from an 
analysis of the inner structure of the echo Tylenda ([2004 ) 
estimated that V838 Mon is at a distance of 8.0 ±2.0 kpc. 
More recently, Crause et al. I|20051 made an analysis of the 
echo expansion observed at the SAAO between May 2002 
and December 2004 and conclude that the distance is 

- 9 kpc. 

Wisniewski et al. ( 20Q3J) have analysed radial velocity 
structures in the interstellar Nal D lines observed in the 
spectrum of V838 Mon. From a comparison with a velocity 
contour map in the Galaxy they conclude that the lower 
limit on the distance is ~ 2.5 kpc. From observations of 
the same lines Kipper et al. ( 2004 ) determined a reddening 
distance > 3 kpc and a kinematic distance > 4 kpc. 

The identification of the B-type companion in the spec- 
trum of V838 Mon in Munari et al. (|2002a|) enables the ap- 
plication of the spectroscopic parallax method. Adopting 
the photomatric data, spectral classification and Eb~v 
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from Munari et al (2002a) as well as the absolute magni- 
tude My — —1.6 for the spectral class B3 V from Schmidt- 
Kaler (1982) one gets a distance to the object of 9.4 kpc. 
(Munari et al., !2UU2al give d = 10.5 kpc, presumably 
adopting a somewhat different calibration of My.) A lower 
limit to the distance from the above data can be obtained 
adopting the ZAMS calibration of Schmidt-Kaler (|1982[) 
which gives My = —1.1 for the B3V type. The distance 
then decreases to 7.5 kpc. From the photometric data of 
Crause et al. I|2()()5|l . Tylenda et al. I|2005b|) derived pa- 
rameters for the B-type companion slightly different from 
those of Munari et al. (|2002a[l . namely V = 16.26 ± 0.02, 
Eb-v = 0.71 and a B4V spectral type. In this case one 
obtaines a distance of 12.4 kpc and 9.4 kpc using the main 
sequence and ZAMS calibrations of Schmidt-Kaler (1982), 
respectively. 

The conclusion of Tylenda et al. (2005b) that 
V838 Mon is very likely to be a young binary system al- 
lows us to assume that the system follows more or less the 
Galactic disc rotation and to estimate the distance from 
the radial velocity of the system. Kolev et al. ( 2002 ) from 
the fairly symmetric emission Ha profile in January 2002 
have derived a heliocentric wavelength of the line of 
6564.155 A which results in a heliocentric radial velocity of 
+62 km s -1 . Kipper et al. pf0l| obtained +59±6 km s" 1 
from emission line profiles in February 2002. The line pro- 
files were rather complex so this estimate is uncertain. 
T. Tomov (2004, private communication) derived a ra- 
dial velocity of ~ +64 km s^ 1 from the H/3 absorption 
line in the spectrum of the B-type component. The above 
values agree very well with the radial velocity of a more 
redshifted component in the interstellar Nal D lines, i.e. 
+64 ± 1 km s" 1 (Zwitter & Munari I2HU21 Kolev et al. 
2002, Kipper et al. 2004). This supports the conclusion 
of Tylenda et al. I|2005b|l that V838 Mon is embedded 
in the interstellar medium. Thus we can adopt that the 
heliocentric radial velocity of V838 Mon is +64 km s _1 . 
Using the results of Dehnen & Binney (199SJ) this can be 
transformed to Vlsr — +53 km s _1 which, adopting the 
Galactic rotation curve of Brand & Blitz (1993|, gives a 
distance of ~ 6.7 kpc. 

Munari et al. ( 2005) from an analysis of interstellar ex- 
tinction, galactic kinematics and properties of the B-type 
companion have concluded that the distance to V838 Mon 
is ~ 10 kpc. 

Summarizing the above distance determinations we 
can conclude that observations of V838 Mon itself and of 
its B-type companion lead to consistent results and that 
the system is most likely at a distance between 5 and 
12 kpc. For the purpose of the present study we assume 
that the distance to V838 Mon is 8.0 kpc. 

From interstellar components of K I and Na I lines 
in the spectrum of V838 Mon Zwitter & Munari l|2002[l 
estimated an interstellar extinction Eb-v — 0.80 ± 0.05 
and concluded that this value would imply a distance ^ 
3 kpc. Munari et al. (2002c|) have, however, argued that 
this value of extinction was too large as their distance 
estimate gave 790 pc and they adopted Eb-v = 0.50. 



This distance estimate was however wrong, as discussed 
above, and V838 Mon is certainly more distant than 3 kpc. 
Therefore the E B -v estimate in Zwitter & Munari (2002) 
can be considered as reasonable. 

Kimeswenger et al. CDD2J have compared effective 
temperatures estimated from spectroscopy with photo- 
metric colours and found that 0.5^i?B-v^0.9. For their 
final discussion Kimeswenger et al. adopted Eb-v — 0.7. 

Kipper et al. (2004) have applied different methods 
for estimating the reddening. From equivalent widths of 
the Na I interstellar lines they derived Eb-v = 0.49. 
Using equivalent widths of diffuse interstellar bands they 
obtained Eb-v between 0.41 and 2.1. From dust distri- 
bution maps these authors estimated an upper limit of 
E B -v ^ 0.89 in the direction of V838 Mon. 

Probably the most reliable estimate of the reddening 
done in Kipper et al. is that resulting from their Table 3. 
These authors have classified their spectra of V838 Mon 
obtained between 4 February and 2 April 2002 using com- 
parison spectra of supergiants. The (B — V) values ex- 
pected from the classification are then compared with the 
observed values giving Eb-v- From our compilations we 
have supplemented the missing data in their Table 3 for 
26 March and 2 April and then from all the results we 
have obtained E B -v = 1-09 ± 0.17. 

As discussed in Tylenda et al. (2005b) the available 
photometric data on the B-type companion give Eb-v 
between 0.7 and 1.0. 

In a recent paper Munari et al. ( 2005 ) have determined 
the interstellar extinction from interstellar Nal and KI 
lines, distribution of the extinction along the line of sight 
and properties of the B-type companion. They conclude 
that E B -v ^ 0.87. 

In the present study we assume Eb-v — 0.9. This 
value is a good compromise between the various results 
obtained for V838 Mon itself and its B-type companion, 
especially if estimates as low as 0.5 are not considered as 
being too low for a distance > 5 kpc near the Galactic 
plane. We also note that the best fits between the photo- 
metric results analysed in the present study and the refer- 
ence supergiant spectra are in most cases obtained when 
Eb-v — 0.9 is adopted. 

4. Results 

Throughout the present paper the epoch of observations 
is given in days counted since 1 January 2002. We devide 
the time span discussed into three phases (for the light 
curve of V838 Mon see e.g. Munari et al. I2002bl 12002^1 
or Crause et al. f2~UU3l IMJS) . The first one, called pre- 
emption, lasted untill 1 February 2002, i.e. it is includes 
epochs < 32 days. The object then slowly evolved between 
V = 10 — 11. The second one, called eruption, lasted untill 
mid-April 2002 covering the epochs between 32 and ~ 
100 days. The object was then brighter than 9 magnitudes, 
twice reaching V ~ 7. In mid- April V838 Mon started a 
fast decline in V, the period with epochs J; 100 days is 
called decline. 
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Fig. 1. Examples of the fits of reference spectra to the observed magnitudes of V838 Mon. Left panel: pre-eruption 
and eruption phases. Right panel: decline phase. Symbols - observed magnitudes. Full curves - supergiant spectra. 
Dotted curves - blackbody distributions. All reference spectra have been reddened with Eb-v = 0.9. The curves are 
labelled with the epoch of observations (days since 1 Jan. 2002). Note that the curves and the data corresponding 
to days 40.8, 47.6 and 72.8 (left panel) have been shifted by —2, —4 and —6 magnitudes, respectively, while those 
corresponding to days 301, 503, 861 and 1075 (right panel) have been shifted by 2, 5, 8 and 10 magnitudes, respectively. 
Parameters of the fits can be found in Table 



Figure Q presents examples of the reference spectra 
(full curves: supergiants - dotted curves: blackbodies) fit- 
ted to the photometric data (symbols). The individual 
spectra have been labelled with the epoch of observations. 

Tabled shows the results from fitting the photometric 
data. Column (1) shows the epoch of observations (days 
since 1 January 2002). Column (2) gives the spectral range 
used in the fitting procedure (e.g. BI stands for BVRI 
while BL means BVRIJHKL) and the references to the 
data (explained in the bottom of Table . The spectral 
type, the corresponding effective temperature and the an- 
gular photosphcric radius (9 is in radians) resulting from 
fitting the standard supergiant spectra to the data are 
presented in columns (3), (4) and (5), respectively. The 
(linear) photosphcric radius and the luminosity, given in 
columns (6) and (7), have been calculated assuming a dis- 
tance of 8 kpc. The second line, labelled BB in column 
(3), for epochs greater than 110 days gives the results of 
the fits using blackbodies as reference spectra. 

The left panel of Fig. ^ displays the spectra during the 
pre-eruption and eruption phases. The U magnitude has 
not been taken into account in the fitting procedure. This 
band is dominated by the Balmcr continuum which is very 
sensitive to nonstandard phenomena such as departures 
from hydrostatic equilibrium, non-LTE and winds. As can 
be seen from Fig. ^ (day 12.2), in pre-eruption the U mag- 
nitude was below that expected for standard supergiants. 



However, during the sharp maximum at the begining of 
February the Balmer continuum was well above that in 
supergiants (see day 35.8 in Fig.^). This period was also 
marked by the appearance of strong Balmer emission lines 
(see e.g. Osiwala 12QQ3J. Also the quality of the obtained 
spectral fit was then less good in the sense that the ob- 
served spectrum was somewhat flatter than the standard 
one. In later epochs, after the sharp maximum, the over- 
all spectral fit became better, including the U magnitude 
(see days 40.8 and 47.6 in Fig. QJ. Thus the object then 
resembled a standard supergiant quite closely. 

As can be seen from Tabled for a number of epochs in 
pre-eruption and eruption the fitting has been done using 
only the BVRI magnitudes, as no JHKL measurements 
were available for these dates. The results from these fits 
are expected to be less reliable. However, as discussed in 
Sect. |S1 there is no significant systematic difference be- 
tween them and those based on the BVRIJHKL data. 

The right panel of Fig. ^ shows examples of the fits 
obtained in the decline phase. As explained in Sect. [21 the 
fits for this period were obtained using both the extrapo- 
lated standard supergiant spectra and the blackbody dis- 
tributions. In all fits in the decline period a contribution 
from the B-type companion, the same as in Tylenda et al. 
(2005b), has been taken into account. This contribution is 
important only in the U BV bands. As can be seen from 
the right panel of Fig.^ the observed spectra change their 
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Table 1. Results from fitting the supergiant and black- 
body spectra (the latter only for day >110) to the photo- 
metric data. See text for explanations of the columns. 
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_B/-c3 


G8.4 


4560. 


8.46 


1230. 


5.9E5 


57.1 


BJ-m2 


K0.6 


4360. 


8.42 


1360. 


6.0E5 


58.1 


BJ-m2 


K0.6 


4360. 


8.41 


1380. 


6.2E5 


58.8 


B/-k2 


G9.3 


4480. 


8.44 


1290. 


6.0E5 


59.8 


BI-c3 


G6.3 


4740. 


8.47 


1200. 


6.6E5 


60.8 


Bl-ci 


G4.4 


4920. 


8.49 


1140. 


6.9E5 


62.2 


BI-m2 


G3.6 


5000. 


8.48 


1170. 


7.7E5 


62.9 


BI-k2 


G3.4 


5020. 


8.48 


1180. 


8.0E5 


63.2 


BI-m2 


G2.1 


5180. 


8.50 


1120. 


8.2E5 



Table Q (continued) 



day 


data 


Sp. 


T e f£ 


-Ig0 


R/Rq 


L/L Q 


63.8 


BI-k2 


G3.0 


5060. 


8.47 


1210. 


8.6E5 


65.8 


BI-k2 


G2.5 


5130. 


8.45 


1260. 


9.9E5 


66.8 


BI-k2 


G2.5 


5130. 


8.44 


1290. 


1.05E6 


67.1 


BI-m2 


G2.7 


5100. 


8.42 


1340. 


1.10E6 


67.8 


BI-k2 


G2.6 


5110. 


8.43 


1320. 


1.07E6 


68.2 


BI-m2 


G2.7 


5100. 


8.42 


1340. 


1.10E6 


68.8 


BI-k2 


G2.9 


5070. 


8.42 


1350. 


1.08E6 


69.1 


BI-m2 


G2.4 


5140. 


8.42 


1330. 


1.12E6 


69.8 


BI-k2 


G3.4 


5010. 


8.41 


1370. 


1.07E6 


71.8 


BL-c3,k2 


G5.3 


4830. 


8.38 


1460. 


1.05E6 


72.8 


BL-c3,k2 


G6.6 


4710. 


8.37 


1520. 


1.03E6 


73.8 


BI-k2 


G8.5 


4560. 


8.33 


1650. 


1.06E6 


85.8 


BI-k2 


K3.0 


4090. 


8.23 


2080. 


1.09E6 


86.8 


BI-k2 


K3.2 


4060. 


8.22 


2120. 


1.10E6 


88.8 


BI-k2 


K3.7 


4000. 


8.21 


2200. 


1.11E6 


90.2 


BI-m2 


K4.1 


3950. 


8.18 


2350. 


1.21E6 


91.8 


BI-k2 


K4.3 


3920. 


8.20 


2250. 


1.07E6 


93.8 


BI-k2 


K4.4 


3910. 


8.20 


2250. 


1.07E6 


98.1 


BI-m2 


K9.8 


3640. 


8.14 


2580. 


1.05E6 


99.8 


BI-k2 


K9.1 


3670. 


8.18 


2340. 


8.9E5 


102.8 


BI-k2 


M1.0 


3550. 


8.20 


2250. 


7.2E5 


119.8 


BL-c3 


M6.7 


2490. 


8.10 


2800. 


2.7E5 




IL 


BB 


3060. 


8.30 


1760. 


2.5E5 


120.7 


BL-cS 


M7.1 


2420. 


8.08 


2920. 


2.6E5 




IL 


BB 


3000. 


8.30 


1790. 


2.3E5 


153.2 


BL-c3 


M9.4 


2020. 


8.08 


2970. 


1.3E5 




IL 


BB 


2070. 


8.16 


2430. 


9.8E4 


232.0 


BL-hm 


Mil 


1780. 


8.05 


3190. 


9.2E4 




IL 


BB 


1710. 


8.06 


3120. 


7.5E4 


254.1 


BL-co 


M9.7 


1980. 


8.16 


2450. 


8.3E4 




IL 


BB 


1840. 


8.14 


2540. 


6.7E4 


301.0 


BK-co 


M8.5 


2180. 


8.26 


1960. 


7.8E4 




IK 


BB 


2120. 


8.27 


1890. 


6.5E4 


385.9 


BK-co 


M8.2 


2230. 


8.32 


1690. 


6.4E4 




IK 


BB 


2210. 


8.36 


1560. 


5.2E4 


391.8 


BK-co 


M8.1 


2240. 


8.33 


1670. 


6.3E4 




IK 


BB 


2220. 


8.37 


1530. 


5.1E4 


474.0 


BK-co 


M7.6 


2330. 


8.40 


1410. 


5.3E4 




IK 


BB 


2390. 


8.47 


1210. 


4.3E4 


503.0 


BK-co 


M7.5 


2340. 


8.41 


1370. 


5.0E4 




IK 


BB 


2380. 


8.47 


1190. 


4.1E4 


619.1 


BK-co 


M6.8 


2470. 


8.48 


1180. 


4.7E4 




IK 


BB 


2580. 


8.55 


1000. 


4.0E4 


650.0 


BK-co 


M6.8 


2480. 


8.49 


1140. 


4.4E4 




IK 


BB 


2520. 


8.54 


1020. 


3.8E4 


861.0 


BK-co 


M6.6 


2510. 


8.53 


1050. 


3.9E4 




IK 


BB 


2630. 


8.60 


890. 


3.4E4 


1075.0 


BK-co 


M5.9 


2620. 


8.64 


800. 


2.8E4 




IK 


BB 


2280. 


8.54 


1020. 


2.6E4 



Data references: c3 - Crause et al. (2003J, c5 - Crause et al. 
( 2005 ), hm - Henden & Munari J2002H . k2 - Kimeswenger 
et al. (1277021 . m2 - Munari et al. (2002cjl 



slope below the I band, i.e. for shorter wavelength bands 
the spectrum is steeper than for the longer ones. This ef- 
fect, due to strong molecular bands significantly lowering 
the flux at shorter wavelengths, is roughly reproduced in 
the supergiant spectra. The blackbody spectra do not take 



this into account. Therefore, when fitting blackbodies we 
excluded the BVR magnitudes from the fitting procedure. 
For the epochs later than 300 days the L magnitude was 
also omitted from the fitting (both using supergiant as 
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well as blackbody spectra) . This band was likely to be af- 
fected by the infrared excess discussed in Sect.EJ Inclusion 
or exclusion of the L magnitude in the fitting procedure 
does not significantly change the fit results. 

As can be seen from Fig. Q] given the crudeness of the 
reference spectra used in the decline phase (extrapolated 
supergiants or blackbodies), the obtained fits to the ob- 
servations are reasonable (although not as good as those 
in pre-eruption and eruption). The values of the fitting 
parameters (T ff and 9 in Table Q obtained using the su- 
pergiant spectra are not significantly different from those 
using the blackbodies in most cases. Only for two spec- 
tra obtained near day 120 does fitting a blackbody give 
a T e g larger by almost 600 K and 9 smaller by a factor 
of ~ 1.6 compared to the results of fitting a supergiant. 
As discussed in Sect. 15.21 the object was probably then 
surrounded by an extended shell in a transition from op- 
tically thick to optically thin. In the case of the last spec- 
trum (day 1075) the blackbody fitting also gave T e g about 
350 K lower and, consequently, a somewhat larger effective 
radius than the supergiant fitting. 

Note that for the reasons discussed above we do not 
consider fits based only on BVRI in the decline phase. 

Lane et al. f20D21 attempted to determine the angular 
size of V838 Mon from interferomatric observations done 
in the K band in November-December 2004. Their result 
of an angular diameter of 1.8 mas is slightly larger than 
our result obtained on day 1075 (~ 1.1 mas). One possible 
reason for this difference is that the uniform disc used in 
Lane et al. to deconvolve the data might be not a good 
model for the photosphere of V838 Mon. Another possi- 
ble reason is that although the K band, according to our 
analysis, is dominated by the central object, a certain con- 
tribution from an extended infrared emission, discussed in 
Sect. El may also be important. This might explain why 
an elliptical model better fits the data of Lane et al. than 
a circular one. It thus would be desirable to make similar 
interferometric measurements but in shorter wavelengths, 
say in / or J, where the central object certainly dominates 
the flux, as well as in longer wavelengths, say in M or N, 
where the expanding matter is dominant. 



5. Evolution of the object 

Figure [3 shows the evolution of the principal parameters 
of V838 Mon with time. The effective temperature, lu- 
minosity and effective radius from Table ^ are displayed 
in the upper, middle and bottom panel, respectively. On 
the abscissa in each panel the epoch of observations (days 
since 1 Jan. 2002), in the logarithmic scale is given. Star- 
like symbols show the results obtained from fitting the 
supergiant spectra to the BVRIJHKL (or BVRIJHK 
for epochs > 300 days) magnitudes, i.e. those labelled with 
BL or BK in column (2) of Tablefl] The same but derived 
from the BVRI magnitudes (labelled BI in column (2) of 
Table Q is represented with full dots. Open symbols dis- 
play the results obtained from fitting the blackbody spec- 
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Fig. 2. Evolution of the effective temperature (upper), 
luminosity (middle) and radius (bottom) of V838 Mon. 
The luminosity and radius are in solar units. The ab- 
scissa gives the logarithm of time in days counted since 
1 Jan. 2002. See text for explanation of the symbols and 
curves (in the middle and bottom panels). 
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tra, i.e. those labelled BB in column (3) of Tabled Dotted 
and dashed curves are discussed Sect. 15.21 and 15.31 

As can be seen from Fig. |3 in pre-eruption and erup- 
tion (log time 2.0), there is no significant system- 
atic difference between the results obtained from fitting 
only the BVRI magnitudes and those derived form the 
BVRIJHKL measurements. The full dots consistently 
follow the same evolutionary trend as the asterisks. In the 
decline phase, given the uncertainties due to the reference 
spectra, the results obtained from fitting the blackbodies 
(open circles) are fairly consistent with those based on the 
supergiant spectra (asterisks). 

In various estimates done later in this paper, requir- 
ing parameters of the central star, we adopt, following 
Tylenda et al. l|2005b|l . standard values for a B3V star, 
i.e. M* = 8 Mq and i?* = 5 Rq. 

5.1. Pre-eruption: epochs earlier than 32 days 

In the pre-eruption phase (log time Ss 1.5 in Fig. |5J) 
V838 Mon seemed to have reached an equilibrium or 
relaxation configuration. This phase followed an initial 
event which had presumably taken place a few days be- 
fore 1 January 2002 (e.g. Munari et al. l2002c|) . The pho- 
tosphcric radius between days 10-32 remained at an al- 
most constant value of ~ 350 Rq, suggesting that the 
photosphere was in a quasi-hydrostatic equilibrium. The 
effective temperature reached a maximum at days 12-15 
and subsequently slowly declined. The luminosity showed 
a similar trend and at the end of the pre-eruption phase 
the objects was half as luminous as at the beginning. 

This behaviour suggests that the initial event, which 
inflated the stellar envelope to ~ 350 Rq, had proper- 
ties of a short-lived energy burst, after which the energy 
sources dropped significantly, provoking a cooling of the 
photospheric regions and a decline in the apparent lu- 
minosity. In principle, cooling and a luminosity decline 
should be followed by a contraction of the photosphere. 
However, for an 8 Mq star and the observed photospheric 
radius the dynamical time is ~ 30 days. Thus the inflated 
stellar envelope did not have enough time to react during 
the pre-eruption phase. 

Within the above hypothesis one can estimate the 
mass of the inflated envelope assuming that during the 
pre-eruption phase it radiated away a significant part 
of its internal energy. The latter can be estimated from 
Eq. 1|A.16J| . Integrating the observed luminosity over the 
pre-eruption phase and assuming (unperturbed) stellar pa- 
rameters, M* = 8 Mq and i?* = 5 Rq, one gets a mass 
of the inflated envelope, M c ~ 4 x 10~ 3 Mq. This value 
should be considered with caution. On the one hand, the 
observed luminosity and temperature evolved only slightly 
during pre-eruption, so only a part of the envelope energy 
radiated away. On the other hand, it is likely the energy 
sources that had presumably initiated pre-eruption did not 
extinguish completely and were contributing to the ob- 
served luminosity in the course of the pre-eruption phase. 



Spectroscopic observations obtained in the pre- 
eruption phase showed numerous lines with P-Cygni pro- 
files indicating significant mass loss. The typical outflow 
velocity was of 200-300 km s _1 with a terminal one reach- 
ing 500 km s" 1 (Munari et al. 20023 Osiwala et al-EDOH- 
This mass loss, observed during the whole pre-eruption 
phase, is not surprising, as the stellar luminosity was then 
of a quarter of the Eddington limit for an 8 Mq star. 

5.2. Eruption: epochs between 33 and ~ 100 days 

The gentle relaxation in the pre-eruption phase was sud- 
denly interupted on 2 February 2002 (day 33). The event 
was very dramatic. For almost 3 days V838 Mon bright- 
ened by a factor of ~ 25 in luminosity while the effective 
temperature increased from ~ 4600 K to ~ 7200 K. The 
event was also followed by a fast rise of the effective ra- 
dius. A major qualitative difference between the evolution 
of V838 Mon in the whole eruption phase compared to 
that in pre-eruption is a more or less steady increase of 
the effective radius in eruption versus a roughly constant 
radius in pre-eruption. At the end of the eruption phase 
the effective radius reached values approaching 3000 Rq, 
compared to ~ 350 Rq in pre-eruption. 

However, as can be seen from Fig. (bottom panel), 
the rise of the effective radius was not monotonia One can 
easily distinguish three periods during which the radius 
was steadily increasing, interupted by two short phases of 
contraction. A simplistic interpretation of this behaviour 
might be that the effective radius evolution reflects the 
behaviour of the stellar envelope, namely that each phase 
of the envelope expansion is followed by a period during 
which the envelope is contracting. However, the expansion 
velocities as estimated below are much higher than the lo- 
cal escape velocities. Thus there is no way to stop the 
expanding matter at a certain time and to force it to con- 
tract. A more likely hypothesis is that during the contrac- 
tion phase the hitherto photospheric layers become, due 
to expansion, optically thin and the photosphere starts to 
form in deeper denser regions. In this case the observed 
evolution of the effective radius during eruption would sug- 
gest ejection of three consecutive shells. 

To futher investigate the evolution of the effective ra- 
dius we have plotted three dotted lines in the bottom panel 
of Fig. [3 showing radii increasing with constant velocity. 
The parameters of the lines, i.e. velocity and epoch of the 
zero radius, are 800 km s _1 and 28.4 days, 400 km s _1 
and 30.5 days, 270 km s _1 and 28.0 days for lines (a), (b) 
and (c), respectively. The fact that the three lines start at 
almost (within 2.5 days) the same time moment strongly 
suggests that the whole eruption phase resulted from a 
single outburst event which took place in the last days of 
January 2002, deep in the stellar envelope already inflated 
by the pre-outburst event. 

The above suggestion that during the eruption phase 
V838 Mon ejected three consecutive shells poses some 
problems when spectroscopic data are considered. The 
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velocity of the first shell (line (a) in the bottom panel 
of Fig. [2J) of 800 km s -1 is much higher than the ob- 
served values. In general, the observed velocities dur- 
ing eruption were similar or even lower than during pre- 
emption. The terminal velocity derived by Kipper et al. 
PMjl on 4/5 Feb. 2002 (days 35-36) was between 240 
and 390 km s _1 . Even if the systematic velocity of the 
object, ~ 65 km s _1 as discussed in Sect. |3 is taken into 
account, this is much lower than the velocity of our hypo- 
thetical first shell. Typical velocities of maximum absorp- 
tion in the P-Cygni profiles in eruption were between 60 
and 300 km s _1 with a general tendecy of decreasing (or 
the appearance of lower velocity components) with time 
(Kolev et al. l2l¥)2l Crause et al. l271T)3l Wisniewski et al. 
2003, Kipper et al. [2004). If the systematic velocity of 
the object is added, these values are consistent with the 
velocities of the shells (b) and (c) in Fig. [21 

A varying photosphcric radius need not reflect a mou- 
vement of the matter in the photosphcric layers. As dis- 
cussed above, the shrinkage of the photosphere, observed 
on days ~ 38 and ~ 60, cannot be considered as due to a 
collapse of the envelope. Similarily an increase of the pho- 
tosphcric radius does not necessarily require expansion of 
a dense, optically thick shell. In the case of an optically 
thick wind, e.g. as discussed in Bath & Shaviv (1976 ) for 
classical novae, the photosphere is formed in the wind and 
its position is mainly determined by the mass loss rate, M. 
If M is increasing the photosphere moves outward, while a 
decrease in M leads to photospheric contraction. However, 
even if M remains constant, but the temperature of the 
wind matter varies, the position of the photosphere would 
also vary because of the varying opacity. We argue that 
this is what happened in the luminosity peak observed 
during the first ~ 10 days of the eruption phase. 

Let us assume a steady-state spherically symmetric 
wind with a mass loss rate, M, and expanding with a 
velocity, v w . Then the density, p, of the wind matter is 
distributed with a radius, r, as 



M 



4-7T r 2 v„ 



(1) 



The effective photosphere is formed at a radius, Rq, sat- 
isfying the condition 



k pair 



Ro 



3' 



(2) 



where k is the opacity. Assuming a constant mass loss 
rate, velocity and opacity above the photosphere, Eq. 
can be integrated, using Eq. (JJJ, leading to 



KpoRo = 



(3) 



where po — M / (4ir R 2 , v w ) is the density in the photo- 
sphcric region. For a given chemical composition, the opac- 
ity, k, is a function of temperature and density. Therefore, 
if values for the effective temperature and radius are given, 
the density, po, can be derived from Eq. (0). Next, using 



4 s - 1 



o -2 




50 100 150 

Lime (days since 1 Jan. 2002) 

Fig. 3. Mass loss rate derived assuming an optically thick 
wind model (see text for more details) - large symbols. 
Small symbols - effective radius (the same as in Fig. [5] 
but plotted as 3 \og(R/R Q ) - 10.5). 



the wind velocity, v w , the mass loss rate M can be calcu- 
lated from Eq. QJ. 

Figure [3] shows the mass loss rate derived from the 
values of the effective temperature and radius listed in 
Tabled using the above steady-state wind model. As the 
opacity we have used the Rosseland means as given by 
Alexander & Ferguson fl9"91| for the chemical composi- 
tion X = 0.7 and Z = 0.02. A constant wind velocity, 
v w = 300 km s -1 , has been assumed. To make the discus- 
sion easier the effective radius has also been plotted in the 
figure with small symbols. 

The values of M plotted in Fig. [3] are reasonable es- 
timates of the mass loss rate only if the wind is optically 
thick (so the effective photosphere is formed in the wind) 



and builds up an 



density distribution above the 



photosphere. The latter condition takes place if the wind 
does not vary significantly on a time scale ~ Rq/v w . This 
time scale is ~ 10 days in pre-eruption. Therefore the very 
fast drop in mass loss rate, by more than 2 orders of mag- 
nitude on a time scale of ^ 1 day, obtained on day 32 is 
not real. If the wind in pre-eruption was optically thick 
and the drop in mass loss rate was real, the density dis- 
tribution above the photosphere would have significantly 
changed after ~ 10 days. However, this was the beginning 
of the eruption phase and the rise in luminosity should 
be followed by an increase in mass loss rate rather than a 
drop. However what is striking in Fig.[3]is a constant mass 
loss rate obtained during the first luminosity peak in erup- 
tion, i.e. between days 32 and 42. During these 10 days all 
the stellar parameters, i.e. luminosity, effective tempera- 
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ture and radius, varied by significant factors. For instance, 
the luminosity first increased by a factor of 25 and then 
dropped by a factor of 2. Yet the mass loss rate was con- 
stant within 10-15%. We therefore suggest that this was 
the time period during which the wind was optically thick 
and our model of a steady-state wind gives a good in- 
terpretation of the evolution of the object. However, the 
infered mass loss rate charactarises the wind during the 
pre- eruption phase rather than in eruption. Thus we in- 
terepret the results obtained in pre- eruption and eruption 
as follows. 

In the pre-eruption phase, as discussed in Sect. l5.il the 
photosphere was formed in layers approximately in hy- 
drostatic equilibrium. The object had a wind as indicated 
by the observed P-Cygni profiles. However the mass loss 
rate was much lower than that shown in Fig. [3] during this 
phase. Instead we argue that M was then close to the value 
derived from the analysis of the first 10 days in eruption, 
i.e. it was ~ 5 x 10 _4 MQ/yr. Note that, as discussed in 
Sect. 15 .J the luminosity of the object in pre-eruption was 
of a quarter of the Eddington limit so a significant mass 
loss could be expected. At the above mass loss rate and a 
photospheric temperature of ~ 5000 K, the wind, during 
the pre-eruption phase, was optically thin and seen only 
in the lines. This roughly steady-state wind built up an 
~ r~ 2 density distribution above the photosphere. On day 
32 a strong luminosity wave (or a series of strong shock 
waves) , presumably generated at the base of the inflated 
envelope, reached the photosphere and quickly heated up 
the regions above it. This resulted in an abrupt increase of 
the opacity (by 2 orders of magnitude if the temperature 
increased from 5000 K to 7000 K - Alexander & Ferguson 
1994) and the hitherto optically thin regions built by the 
wind became optically thick. This was observed as a very 
fast expansion of the photospheric radius leading to an 
apparent velocity of ~ 800 km s" 1 , much faster than the 
wind velocity observed in the lines. Subsequent cooling of 
the wind matter decreased the opacity, which resulted in a 
contraction of the effective photosphere observed between 
days 36 and 42. 

The luminosity wave was followed by an intense mass 
loss (the observed luminosity between days 33 and ^110 
was above the Eddington limit for an 8 Mq star) which 
reached the shrinking photosphere near day 42, as indi- 
cated by a large increase of M in Fig. |21 Note, however, 
that the obtained values of M should again be regarded 
with caution. The assumption of an r~ 2 density distribu- 
tion in the model wind is not satisfied if the mass loss rate 
is rapidly changing, as seen in Fig. |3| after day 42. In this 
case the obtained evolution of M in the later epochs, when 
combined with the evolution of the effective radius, seems 
to suggests that mass loss following the fast luminosity 
wave took place in two shells. These shells were presum- 
ably formed by the initial eruption event which took place 
at the end of January, as discussed above, and followed 
the initial strong luminosity wave. Interaction of the lu- 
minosity wave with the envelope inflated in pre-eruption 



also might have been important for the formation of the 
shells. 

The first shell, somewhat faster than the second one, 
became visible near day 42, when it met the photosphere 
contracting in the pre-eruption wind. Subsequently, it was 
observed as an expanding and cooling photosphere radi- 
ating at a roughly constant luminosity of ~ 6 x 10 5 Lq. 
Due to expansion and cooling this shell became partly 
transparent around day 58 and the second shell started 
to be visible. This resulted in a shrinkage of the effec- 
tive photosphere observed between days 58 and 63. At 
the same time the photosphere became hotter while the 
luminosity started to increase and quickly reached a value 
of ~ 1 x 10 6 Lq remaining more or less constant between 
days 66 and 100. Judging from the evolution of the effec- 
tive radius and luminosity, the second shell became partly 
transparent between days 100-120. The discrepancy in the 
effective radius derived from fitting supergiant spectra and 
blackbodies on day 120 (see Tableland Fig.[2J) can be ex- 
plained by the presence of an extended, partly transparent 
envelope. At a temperature of 2500-3000 K the matter is 
opaque in the optical due to molecular absorption so the 
effective radius derived from the BVR bands (important 
when fitting a supergiant spectrum) would correspond to 
the outer radius of the envelope. In the infrared the enve- 
lope is more transparent and it is possible to see deeper 
and hotter layers. Thus a higher effective temperature and 
a smaller effective radius is obtained from fitting a black- 
body to the IJHKL measurements. 

From the parameters observed when a shell becomes 
transparent one can estimate the mass of the shell, as the 
optical thickness of the shell is then expected to be 

2 

t = k P AR s ~- (4) 

where Ai? s is the radial thickness of the shell. If AR S , 
temperature, T, and k(T, p) are known, the density, p, can 
be derived from Eq. ijljl. Then, knowing the shell radius, 
R s , the mass of the shell, M s , can be obtained. 

As discussed above the first shell started to become 
partely transparent on day ~ 58. We can thus assume the 
effective radius and temperature obtained at this epoch 
(see Table ^| as estimates of the shell radius and tem- 
perature, i.e. R s = 1.3 x 10 3 Rq and T = 4500 K. 
From a difference between curves (b) and (c) in the bot- 
tom panel of Fig. |2 at this epoch, we can estimate that 
AR S ~ 0.2 R s . Adopting the Rosseland mean opacity as 
tabulated in Alexander & Ferguson Q1994[l we finally ob- 
tained M s ~ 0.08 Mq. 

The second shell presumably became transparent be- 
tween days 100-120. Assuming that Eq. Q was satisfied 
on day 100 we have (see Table [Q R s = 2.5 x 10 3 Rq and 
T = 3500 K, which results in M s ~ 0.45 Mq. If we assume 
that the shell became partly transparent on day 120 and 
that the results obtained from fitting a supergiant spec- 
trum (see Table nj are estimates of the shell parameters, 
i.e. Rs = 3.0 x 10 3 Rq and T = 2500 K, the result is 
M s ~ 0.65 Mq. 
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The above estimates show that the second shell was 
~ 7 times more massive than the first one. However, 
the obtained absolute values probably overestimate the 
real shell masses. The reason is that the Rosseland mean, 
which very well approximates the true opacity in opti- 
cally thick conditions, is likely to underestimate the real 
opacity in partly transparent media, particularly in con- 
ditions where the Planck mean is much higher than the 
Rosseland mean. This is the case for the temperature and 
density range considered above where, according to the 
results of Alexander & Ferguson l|1994f> . the Planck mean 
is ~ 2 orders of magnitude higher than the Rosseland 
mean. If the Planck mean is used instead of the Rosseland 
mean in Eq. (TJJ the results are M s ~ 8 x 10 -4 Mq and 
~ 4 x 10~ 3 Mq for the masses of the first and second shell, 
respectively. These values can be considered as lower lim- 
its to the shell masses. 

5.3. Decline: epochs later than ~ 100 days 

This phase is first marked by a continuous decline in lu- 
minosity. This decline was initially fast. For about 50 days 
V838 Mon faded by an order of magnitude. Next the time 
scale became longer and after almost three years the ob- 
ject was ~ 40 times less luminous than in eruption. 

The initial decline, untill day ~ 230, was also char- 
acterised by a decreasing effective temperature which 
dropped to ^ 2000 K. The results of our analysis during 
this phase are particularly uncertain. The spectral fitting 
was based on an uncertain extrapolation well beyond the 
range of the standard supergiant spectra. Also the black- 
body spectra were not successful in fitting the optical and 
infrared magnitudes at the same time. The evolution of 
the effective radius was particularly uncertain, as can be 
seen from the bottom panel of Fig. [21 

In the later decline, i.e. days 250, the object shrunk 
in effective radius. Also the trend in the effective temper- 
ature was reversed. The object, although still very cool 
for a star, slowly became hotter. If the observed evolution 
of V838 Mon is plotted on the HR diagram the conclu- 
sion is clear: in the later decline the object declined along 
the Hayashi track. The only known stellar eruptive object 
that shows a similar evolution is V4332 Sgr (Tylenda et al. 
2005a) which is usually considered to be of the same class 
as V838 Mon. 

In the stellar evolution the only case of decline along 
the Hayashi track is the protostellar phase when the main 
source of energy is gravitational contraction. This strongly 
suggests that V838 Mon in the late decline also is powered 
by a gravitational collapse of its envelope inflated during 
eruption. As shown in Tylenda ct al. (2005a) this hypoth- 
esis can satisfactorily explain the observed evolution of 
V4332 Sgr for almost 10 years after its eruption. Within 
this hypothesis one can estimate the mass of the inflated 
envelope. 

Appendix fX] presents an approach to the problem of 
a gravitationally contracting envelope. The structure of 



the envelope is treated in a polytropic approximation 
which seems to be better than the r -5 / 2 density distri- 
bution assumed in Tylenda et al. (2005a). For the effec- 
tive temperatures observed during decline it is reasonable 
to assume that the envelope is convective so n = 3/2. 
The dashed curves in the middle and bottom panels 
of Fig. [21 were derived from a numerical integration of 
Eq. (|A~20)l (using Ecis. lA"T2l and lA~l"3)l . We found that 
T e s = 6120 — 1200 \og(R/ R@) is a good approximation to 
our results obtained from fitting the supergiant spectra to 
the observations for days ^ 250. This relation was used in 
our modelling of the envelope contraction. M* = 8 Mq 
and i?» = 5 Rq have been assumed as the parameters 
of the central star. The mass of the envelope was chosen 
to fit the observed radius and in the case of the dashed 
curves in Fig. H it is M c = 0.22 Mq. 

6. Infrared excess 

Infrared spectra of V838 Mon have been analysed 
by Lynch et al. (2004). The observations taken in 
January 2002 (i.e. in pre-eruption) they have been able to 
fit with a single blackbody of 2700 K. They, however, note 
a possible small excess near 10 /jm. The data obtained a 
year later, between December 2002 and February 2003 (i.e. 
in decline), show a more complex nature. In particular, the 
spectrum shows strong absorption bands due to numerous 
molecules characteristic of an oxygen-rich (C/O < 1) mat- 
ter. Lynch et al. ( 2004) conclude that the overall spectrum 
at short wavelengths can then be fitted with a 2000 K 
blackbody, whereas at longer wavelengths (8-13 /jm) the 
continuum indicates 650-800 K. 

We have collected the available measurements in the 
range of A 5 /im. With the shorter wavelength photom- 
etry (discussed in Sect. 0} this allowed us to investigate 
the presence of an infrared excess for six epochs, as given 
in column (1) of Tabled Two of them are in pre-eruption, 
one in eruption and three in decline. The data and results 
are summarized in Table [21 and Fig. 0] 

Standard supergiant colours are poorly known for the 
MN bands. Therefore in order to determine the stellar 
contribution for A ij 5 /im we have used blackbody dis- 
tributions fitted to shorter wavelengths. For the decline 
period these are simply the same blackbody spectra as 
those discussed in Sect. 0] and whose parameters are given 
in Tabled They are shown with dotted curves in the right 
panel of Fig. 0] . In the case of pre-eruption and eruption 
we have fitted blackbodies to the observed magnitudes 
shortward of 5 /im. They are shown with full curves in the 
left panel of Fig. 01 The values of the stellar blackbody 
temperatures, Tbb,s, are given in column (3) of Tabic 
Note that the blackbody temperatures obtained for the 
first three epochs in Table El are not significantly different 
from the values of T e g obtained from fitting supergiant 
spectra (see Table 0. 

As can be seen from Fig. 0] (left panel) the data in pre- 
eruption and eruption show a 0.5-1.0 magnitude excess at 
~ 10 /im. This is in agreement with the analysis of Lynch 
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Fig. 4. Infrared excess in pre-emption and eruption (left panel) as well as in decline (right panel). Symbols - 
observed magnitudes. Full curves - fitted spectra. Dotted and dashed curves (right panel) - stellar and dust spectral 
components, respectively. The curves are labelled with the epoch of observations (days since 1 Jan. 2002). A downward 
arrow for day 258 shows the upper observational limit at 18 /im. Note that the curve and the data on day 15.0 (left 
panel) have been shifted by 2 magnitudes, while those corresponding to days 395 and 650 (right pannel) have been 
shifted by 4 and 8 magnitudes, respectively. The data sources and the parameters of the fitted spectra can be found 
in Tabled 



Table 2. Infrared excess in V838 Mon: data and results. 
See text for explanations of the columns. 



day 


data 


Tbb.s 


TBB,d 


-Ig0 d Ra/Rq L d /L a 


15. 


BVRI(clk3) 
LMN (lrr4) 


5200. 






28.5 


BVRI(kls2) 

JHK{\xrA) 

9-12/mi(klk2) 


4650. 






42. 


C/-L(mhk2) 
10-12pm(klk2) 


6880. 






258. 


B-L(clm5) 
8-18/xm(wc2) 


1840. 






395. 


B-L(clm5) 
J-NQitA) 


2220. 


800. 


7.75 6.3E3 0.29 


650. 


B-L(clm5) 

LAfiV(lrp3) 

9-19/mi(tp3) 


2520. 


600. 


7.43 13.1E3 0.54 



Data references: clk3 - Crause et al. I|2()l)3jl . clm5 - Crause 
et al. (HOB, klk2 - Kaeufl et al. j2"UU2l . kls2 - Kimeswenger 
et al. lUnOU, lr P 3 " Lynch et al. llS0U3"l . Irr4 - Lynch et al. 
dSUl, mhk2 - Munari et al. (|2()()2c^ . tp3 - Tapia & Persi 
JSUUSll . wc2 - Watson & Costero ll2lK72l 



et al. I)2()()4|l of their January 2002 spectrum. 1 As noted in 



Lynch et al. the 10 /im excess is likely to be due to silicate 
emission. 

A similar conclusion can also be drawn from a decline 
spectrum of day 258 (mid-September 2002). For this epoch 
we have more detailed measurements near 10 /im from 
Watson & Costero |2002). As can be seen from Fig. 0] 
(right panel), the stellar blackbody fits the data well up to 
~ 8 /im. Above this wavelength the observed flux steeply 
increases and between 9-13 /im it stays ~ 0.8 magni- 
tude above the stellar continuum. This feature can in- 
deed be interpreted as due to optically thin emission 
from silicate dust. To show this we have added a compo- 
nent proportional to the dust emission coefficient, i.e. to 
Qabs,i/ -By(Td), where Q a bs,^ is a monochromatic absorp- 
tion coefficient for " astronomical silicates" as tabulated in 
Draine ( 1985}. This component is shown as a dashed curve 
on day 258 in Fig. The full curve presents the sum of 
the stellar blackbody and the dust contribution and, as 
can be seen from the figure, it fits the observations very 
well. The dust temperature, Td, has been assumed to be 
800 K as this is the expected temperature of silicate dust 



1 The difference in the blackbody temperature between 
Lynch et al. (200IJ and the present study (2700 K versus 



~ 5000 K) is simply due to the interstellar reddenning cor- 
rection, i.e. no reddenning correction in Lynch et al, versus 
Eb-v = 0.9 in the present study. Indeed, if we apply our fit- 
ting procedure with Eb-v = 0.0 to the JHKLM magnitudes 
of Lynch et al. in January 2002, the resultant temperature is 
~ 2700 K. 
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condensation. This value is also close to the blackbody 
temperature estimated in Lynch et al. (2004| and in the 
present study from the observations in December 2002 - 
January 2003 (day 395, see below). The value of T d can- 
not be precisely constrained from the observational data 
on day 258, however it cannot be lower than ~ 400 K as 
then the predicted brightness at 18 //m is above the ob- 
servational upper limit (shown by a downward arrow in 
Fig. mi. 

We can therefore conclude that in pre- eruption, erup- 
tion and early decline V838 Mon showed an infrared excess 
at ~ 10 /im. This excess was not significant for the global 
energetics of the object and was likely to be due to silicate 
dust condensation in the wind. The dust was probably op- 
tically thin, suggesting a moderate mass loss rate (at the 
distance of dust condensation). Most likely, following the 
discussion in Sect. 15. 21 this concerns the pre-emption wind 
matter expanding above the two dense shells ejected at the 
onset of eruption. 

In later epochs the situation changed considerably. On 
days 395 (data from observations in December 2002 and 
January 2003) and 650 (October 2003) the infrared ex- 
cess increased and cannot be fitted by an optically thin 
silicate emission. In the former epoch Lynch et al ( 2004 ) 
note a silicate feature at ~ 10 /im but it is superimposed 
on a blackbody of 650-800 K. In the latter epoch, as noted 
in Lynch et al. (|2003|) . the feature became even less pro- 
nounced with a continuum emission dominating the ob- 
served spectrum. 

We have been able to fit the whole spectrum of 
V838 Mon on days 395 and 650 by a sum of two blackbody 
distributions, as shown in Fig. 0] (right panel). The param- 
eters of the fits can be found in Table [21 The cooler com- 
ponent (dashed curves in the figure) can be interpreted 
as due to an optically thick dusty region reprocessing the 
radiation from the central object. The fact that we see 
the central object signifies that dust does not significantly 
obscure it. In other words, the covering factor of the dusty 
region is < 1, in accordance with the discussion below. 

As can be seen from Table the dust temperature, 
?BB.d, decreased while the radius, Rd, increased between 
the two epochs. This can be easily interpreted as due to 
expansion of the dusty region. 

The values of log 9d and Rd given in Table |21 are lower 
limits to the dust distance from the central object, as they 
have been derived assuming a radiating sphere. The values 
of Ld/L s given in the last column of Table El can be con- 
sidered as an estimate of the covering factor of the dusty 
region. If this is the case, the dust distance should be in- 
creased by the inverse square root of the covering factor 
and Rd becomes 1.17 x 10 4 Rq and 1.78 x 10 4 Rq on days 
395 and 650, respectively. These values of Rd and the de- 
rived temperature values closely follow the expected rela- 
tion T BB ,d/ T BB,s = {R s /Rd) 1/2 (R s taken from TableHJ. 
Assuming that, as discussed in Sect. 15.21 the dominant 
mass loss was initiated on day ~ 30 one can estimate a 
mean expansion velocity from the above values of Rd- The 
result is 260 km s _1 and 230 km s _1 for days 395 and 650, 



respectively. We can therefore interpret the strong infrared 
excess observed in the later decline as due to dust formed 
in the second, more massive, shell ejected in eruption, as 
discussed in Sect. 15.21 From the condition that the shell 
should be optically thick we can estimate a lower limit 
to its mass. Following Alexander & Ferguson fT9"g4]|l the 
opacity k ~ 2.7 for T ~ 700 K (the Rosseland mean, but 
the Planck mean is not significantly different). Taking the 
shell radius of 1.8 x 10 4 Rq, a radial thickness of 20% 
and the filling factor of 0.5 the condition r 1 implies 
M s > 2 x 10~ 3 Mq. This value is consistent with the 
estimates obtained in Sect. 15.21 (4 x 10 -3 - 0.6 Mq). 

The high value of T BB ,d in Dec. 2002 - Jan. 2003 
(800 K) suggests that dust was being formed during this 
epoch. Probably the whole process of dust formation was 
completed in a somewhat later epoch. This might be the 
reason why the covering factor as measured by L&/L s on 
day 650 was larger than that on day 395. 

7. Summary 

We present the following description of the eruption of 
V838 Mon observed in 2002. 

The eruption consisted of two distinct phases, pre- 
emption and eruption. Both phases are likely to have re- 
sulted from rather short-lived but strong energy bursts. 

The first one, which produced the pre- eruption, oc- 
cured before the discovery of the object so we have no 
data for it. Presumably it took place in the last days of 
December 2001 and its main effect was to produce an enve- 
lope inflated to ~ 350 Rq. In January 2002 the object was 
in a relatively quiet phase, radiating at ~ 7 — 8 x 10 4 Lq. 
The slowly decreasing effective temperature and luminos- 
ity suggest relaxation after the initial event. The object 
had a wind that built up an optically thin envelope, ex- 
tending to ~ 1500 Rq in last days of January 2002. 

Presumably, another eruption event occured on 28- 
30 January at the base of the inflated envelope. This sec- 
ond energy burst was much stronger than the first one 
and generated the main phase, eruption. As a result an 
energy wave, probably in the form of shocks and/or a ra- 
diation wave, was formed, followed by a massive outflow 
of matter. The energy wave reached the photosphere on 
2 February and heated up the extended envelope built up 
by the pre-eruption wind. This increased the opacity and 
made the inner wind layers optically thick, which was ob- 
served as a rapid expansion of the effective photosphere. 
In 3 days the effective radius increased to ~ 700 Rq, 
while the luminosity reached ~ 1 x 10 6 Lq. Subsequent 
cooling of the wind envelope resulted in a decrease in lu- 
minosity and a contraction of the photosphere. On 10- 
11 February the first shell of mass outflow met the con- 
tractiong photosphere and became visible. This first, less 
massive shell, likely to have been produced from the pre- 
eruption envelope compressed by the second energy burst, 
was observed as an expanding and cooling photosphere 
untill ~ 27 February. Then it became partly transparent, 
so that deeper and hotter layers started to be visible. In 
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this way the main mass outflow was observed between 
4 March and ~ 10 April as an initially hotter, but then 
expanding and cooling, photosphere radiating at a lumi- 
nosity of ~ 1 x 10 6 Lq. Thus we show that the evolution of 
the object in the eruption phase, i.e. in February - mid- 
April, was due to a single event rather than a series of 
outbursts as often suggested in the literature. 

Because of no significant energy supply from internal 
sources the object declined in luminosity when the shell of 
the main outflow became transparent and radiated away 
most of its internal energy. A certain part of the mat- 
ter involved in the eruption did not gain enough energy 
to overcome the gravitational field of the central star. It 
formed an inflated envelope whose gravitational contrac- 
tion was the main source of the luminosity observed in the 
decline phase. We estimate that the mass of this envelope 
is ~ 0.2 Mq. 

At the same time the matter ejected during eruption 
expanded and cooled, so that dust formation was possible. 
As a result a strong infrared excess at A <; 3 /im was 
observed since October-November 2002. We estimate that 
the total mass lost from V838 Mon is < 0.6 Mq. This is 
rather a conservative upper limit. The true value is likely 
to be significantly lower. 

We postpone the discussion of the present results in 
terms of possible eruption mechanisms to another paper. 
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by the Polish State Committee for Scientific Research. 

Appendix A: Polytropic envelope 

Suppose that a star of a mass, M*, and a radius, 
has an envelope of mass, M e , inflated to a radius, R c . 
The envelope is in hydrostatic equilibrium. Suppose also 
that M e <C M*, so the gravitational potential within the 
envelope, V(r), is 



V(r) 



GAh 



(A.l) 



Let us assume a polytropic equation of state, i.e. 

P = Kp 1+ ^/ n \ (A.2) 

where P is pressure, p - density, K - polytropic constant, 
n - polytropic index. (A detailed theory of polytropes can 
be found in e.g. Chandrasekhar IT95 71 ) Then the standard 
equation of hydrostatic equilibrium can be integrated giv- 
ing 



^P GM.« GM* 

(n+ 1)~ = — 

P r R c 



(A.3) 



Eq. (|A.3|I can be used with Eq. I|A.2(I to obtain the density 
distribution, i.e. 



P = Po 



1-Xo 



1-x 



where x = r/R Cl xq = R*/R e , and po = p(xo) is the den- 
sity at the base of the envelope. The mass of the envelope, 
Af e , is 



Re 

R. 



M e = / 4irr 2 pdr = 4tt R% p ' '"' 



1 -Xo 



Im(n,x a ) (A. 5) 



where 



I m (n,x )= / (1 - x) n x 2 ~ n Ax. (A.6) 
The potential energy of the envelope, il, is 

l-Re 

O = / Aur 2 V(r) pdr 
Jr. 

= -AnGM*Rlp (-^-\ I e (n,x ), (A.7) 
V 1 - x J 



where 

I c (n,x ) = [ (1 - x) n a; 1 "™ dx. 

Jxq 

Using Eq. (I A . 5|) . Eq. l|A.7Jl can be rewritten as 
G M* M e J e (n,s ) 



Q = — 



Re I m {n,x )' 



(A. 



(A.9) 



According to the virial theorem the internal energy of the 
envelope, U, is 



u = —n. 

2 

Thus the total energy of the envelope, E, is 
G M* M e 7 e (n,a;o) 



E = n + u 



2R C I m (n,x )' 



(A.10) 



(A.ll) 



In astrophysical applications two values of the poly- 
tropic index, i.e. n = 3/2 and n — 3, are of particular 
interest. The integrals in Eqs. ljA.61) and (|A.8|I can be an- 
alytically evaluated in these cases giving 

/ m (3/2,x ) = ^ - - arcsin(a;o /2 ) - - xl /2 (1 - x Q ) 1/2 

-l^a-^+i^a-xo) 5 / 2 , (A.i2) 

J e (3/2, x ) = 1 7r ~ j arcsin(a;o /2 ) - ^ x^ 2 (1 - x ) 1/2 



(A.13) 



(A.4) e 



11 3 1 

/m(3,x ) = -— - In .T + 3x ~ ^^o + 3^0- ( A - 14 ) 

13 1 
7 e (3, x ) = h - + 3 In x a - 3x + -x . (A.15) 

Xq Z Z 

Note that in the case of extended envelopes (x <C 1) 
the envelope energy becomes 

3 G M* M e 



Rn 



for n = 3/2 



(A.16) 
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and 



E 



GM*M e 
2 i?„ \n(R c /R* 



for 



(A.17) 



Assuming that the luminosity of the envelope, L, is due 
to gravitational contraction of the envelope and neglecting 
mass loss from the envelope, we can write 



L = — 



dE 



dt 

G M* M c 
2R n 



(A.18) 



I e (n,x )\ I c (n,x ) d 



dt \I m (n,x )J I m (n,x ) dt 



(In .Re) 



On the other hand, we have 



cir- 



(A.19) 



where T e g is the effective temperature. Combining Eq. 
I|A.19|) and Eq. (|A.19jl we obtain 



(A.20) 



A(lni? c ) = Afln^4 
dr c > dt \ I m (n,x ) 



inR 3 c aT^ s ( I e (n,x ) 



GM*M e 



7 m (n,x ) 



If the behaviour of T c g is known, Eq. (|A.20jl can be numer- 
ically integrated to obtain the evolution of the envelope 
radius. 
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